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One of the most relevant manifestations of the Beresinskii-Kosterlitz-Thouless transition occurs 
in quasi-two-dimensional superconducting systems. The experimental advances made in the last 
decade in the investigation of superconducting phenomena in low- dimensional correlated electronic 
systems raised new questions on the nature of the BKT transitions in real materials. A general 
issue concerns the possible limitations of theoretical predictions based on the XY model, that was 
studied as a paradigmatic example in the original formulation. Here we review the work we have 
done in revisiting the nature of the BKT transition within the general framework provided by the 
mapping into the sine-Gordon model. While this mapping was already known since long, we recently 
emphasized the advantages on such an approach to account for new variables in the BKT physics. 
One such variable is the energy needed to create the core of the vortex, that is fixed within the 
XY model, while it attains substantially different values in real materials. This has interesting 
observable consequences, especially in the case when additional relevant perturbations are present, 
as a coupling between stacked two-dimensional superconducting layers or a finite magnetic field. 

PACS numbers: 74.20.-z, 74.25. Fy, 74.78.Fk 
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I. INTRODUCTION 



Almost 40 years after the pioneering work by Beresinskii 1 and Kosterlitz and Thouless 2 ' 3 the Beresinskii-Kosterlitz- 
Thouless (BKT) transition is still the subject of an intense experimental and theoretical research. On very general 
grounds, what makes it so fascinating is the possibility of a phase transition that is not driven by the explicit breaking 
of a given symmetry, but is based on the emergence of a finite (and measurable) rigidity of the system. The BKT 
transition was originally formulated within the context of the two-dimensional (2D) Xy-model, which describes the 
exchange interaction between classical two-component spins with fixed length S = 1 : 



H XY = -jJ2cos(6 l -6 J ), (1) 

. (ij) 

(N ' .... ........ 

where J is the spin-spin coupling constant and 0i is the angle that the i-th spin form with a given direction, and i are the 
sites of a square lattice. This model admits a continuous U(l) symmetry (encoded in the transformation 0, — > 0j + x), 
that cannot be broken at finite temperature by an average magnetization (S) = (e 10 ) different from zero because of the 
Mermin- Wagner theorem. Nonetheless, the system can become "stiff" at low temperature with respect to fluctuations 
• • . of the variable, leading to a power-law decay of the spin correlation functions, i.e. (e l ( e W -9 WD) ^ r -T/2u-j ^ j n 
contrast to the exponential one expected in the truly disordered state. Such a change of behavior cannot be "smooth", 
i.e. a phase transition occurs in between, which appears to be controlled by the emergence of vortex-like excitations. 
The original argument used if Ref. 2 to capture the temperature scale of such a transition is rather intuitive: in two 
dimensions both the energy E and the entropy S of a single vortex excitation depend logarithmically on the size L of 
the system, so that the free energy reads: 

F = E-TS= (tt,J-2T) In-. (2) 

a 



As a consequence, at temperatures larger than 



Tbkt ^ — , (3) 



free vortices start to proliferate and destroy the quasi-long-ranged order of the correlation functions. 

As it was observed already in the original papers 2,3 , several physical phenomena are expected to belong to the same 
universality class than the XY model (1), as for example the superfluid transition in two dimensions. Later on it was 
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realized that the same idea can be applied also to superconducting (SC) thin films 4,5 , even though in the charged 
super luid the logarithmic interaction between vortices is screened by the supercurrents at a finite distance A = X 2 /d, 
where A is the penetration depth of the magnetic field and d is the film thickness. 6 In practice, for sufficiently thin 
films with large disorder (so that A is also large) the electromagnetic screening effects are weak enough to expect the 
occurrence of the BKT transition. 

As a matter of fact, the case of SC thin films represented one of the most studied applications of the BKT physics. 
In principle, in this case one has also several possibilities to access experimentally the specific signatures of the BKT 
physics. For example, by approaching the transition from below, the superfluid density n s is expected to go to zero 
discontinuously at the BKT temperature Tbkt, with an "universal" relation between ti s (Tbkt) and Tbkt itself 7-9 , 
that is the equivalent of the relation (3), since J oc n s (see Eq. (5) below). Approaching instead the transition from 
above one has in principle the possibility to identify the BKT transition from the temperature dependence of the 
SC fluctuations. Indeed, in 2D the temperature dependence of both the paraconductivity Act = <t — <t„ and the 
diamagnetism \d is encoded in the SC correlation length £(T), 5 which increases approaching the transitions due to 
the increase of SC fluctuations: 

Within BKT theory £(T) diverges exponentially at Tbkt as ^bkt(T) ~ ae 6 /^, 3 ' 8 where t = T/T B kt - 1, 

in contrast 

to the power- law £gl(T) ~ 1/(T—T C ) expected within Ginzburg-Landau (GL) theory 10 . As a consequence, by direct 
inspection of the paraconductivity or diamagnetism near the transition one could identify the occurrence of vortex 
fluctuations, which lead to an exponential temperature dependence of the SC correlation length. 

Quite interestingly, very direct measurements of the BKT universal jump in the superfluid density of SC films 
became available only recently, 11-16 due to the improvement of the experimental techniques, triggered mostly by 
the investigation of high-temperature superconductors in the late nineties. In particular, the use of the two-coil 
mutual inductance technique 17 turned out to be crucial to obtain the absolute value of the superfluid density at zero 
temperature, which is needed to compare the experimental data with the BKT predictions. Recently a great deal 
of information has come also from Tera- hertz spectroscopy 18-20 , which probes the finite- frequency analogous of the 
superfluid-density jump. At the same time, in the last decade new 2D or quasi-2D SC systems emerged where the 
BKT transition is expected to occur. To this category belong for example the nanometer-thick layers of SC electron 
systems formed at the interface between artificial heterostructures made of insulating oxides as LaAlOa/SrTiOa 21 ' 22 
or LaTiOs/SrTiOa 23 , or at the liquid/solid interface of field-effect transistors made with organic electrolytes 24 . A 
second remarkable example is provided by layered 3D systems as cuprate high-temperature superconductors, where 
the weak interlayer coupling makes it plausible that at least in same regions of the phase diagram a BKT transition 
could be at play 25,26 . Even though the existence of a BKT transition in bulk (3D) samples is still controversial, 
as we shall discuss below, nonetheless in cuprates the proximity of the SC phase to the Mott insulator leads to 
a large penetration depth without need of introducing strong disorder, making in principle thin-films of cuprate 
superconductors the best candidate to study BKT physics 12 ' 13 ' 16 ' 20 ' 27,28 . Recently, much attention has been devoted 
also to artificial heterostructures made of cuprates at different doping level 29,30 , where the observation of BKT physics 
in the SC films is complicated even more by the proximity to the non-SC correlated insulator. 

A common characteristic of the cases mentioned above is that the BKT transition is expected to occur in systems 
where electronic correlations are not necessarily in the weak-coupling limit. This can be due to the presence of 
strong disorder, as it is the case for thin disordered films of conventional superconductors, to the artificial spatial 
confinement, as in the SC interfaces, or to the intrinsic nature of the system, as it occurs in cuprate superconductors. 
As a consequence, several experimental results seem to point towards a kind of "unconventional" BKT physics, which 
needs to be addressed using a wider perspective than the one proposed in the original formulation. A paradigmatic 
example of an apparent failure of the standard BKT approach is posed by recent measurements of the universal 
superfluid-density jump in InO 18 ' 19 and NbN 14 ' 15 films. In a quasi-2D superconductor of thickness d the energy scale 
corresponding to the coupling J of the XY model (I) is the so-called superfluid stiffness, which is connected to the 
(areal) density of superfluid electrons p 2 s d = n s d, which in turn is measured via the inverse penetration depth A of the 
magnetic field: 

tfpf _ h*<?d 

4m 167re 2 A 2 ' 1 ' 

This coupling has itself a bare temperature dependence J(T) due to the presence of quasiparticle excitations: however, 
one would expect that at a temperature scale corresponding to the relation (3) free vortices start to proliferate, so that 
n s (T) jumps discontinuously to zero. In the experiments of Ref. 14,15 one can clearly see that as the film thickness 
decreases n s (T) starts to deviate abruptly from its BCS temperature dependence. However, such a deviation seems 
to occur at a temperature lower than the one predicted by Eq. (3). The same observation holds for finite- frequency 
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measurements of n s (T) in InO films 18,19 , casting some doubt on that "universal" relation between the superfluid 
density and the critical temperature that is one of the hallmarks of the BKT transition. 

It is worth noting that while the measurement of the superfluid-density behavior gives access to the most straight- 
forward manifestation of BKT physics, its identification via SC fluctuations is much more subtle. Indeed, according 
to the general result (4), one needs in this case a controlled procedure to first extract the SC fluctuations contribution, 
and then to fit it with the BKT expression for the correlation length. Such a procedure is in general applied to the 
paraconductivity, even though much care should be used to disentangle GL from BKT fluctuations, as it has been 
discussed in a seminal paper long ago by Halperin and Nelson 5 (HN). Indeed, the BKT fit should be applied only 
in the region between the BCS mean-field temperature T c and the true Tbkt, that do not differ considerably in 
thin films of conventional superconductors. In contrast, recent applications to paraconductivity measurements in thin 
films of cuprate superconductors 28 or in SC interfaces 21 ' 31 seem to suggest that in these systems the whole fluctuation 
regime above Tbkt is dominated by BKT vortex fluctuations, and deviations only occur near the transition because 
of finite-size effect. Also in this case, one would like to distinguish unconventional effects due possibly to the nature 
of the underlying system, from spurious results due to an incorrect application of BKT theory. As we discussed 
recently in Ref. [32] a BKT fit of the paraconductivity must be done taking into account from one side the existence of 
unavoidable constraints on the values of the fitting parameters, and from the other side the existence of inhomogeneity 
on mesoscopic scales, that can partly mask the occurrence of a sharp BKT transition. A very interesting example 
of application of such a procedure has been recently provided by NbN thin films 15 , where the direct comparison 
between superfluid-density data below Tbkt and resistivity data above Tbkt provided a paradigmatic example of 
BKT transition in a real system. 

The issue of the inhomogeneity emerges also within the context of high-temperature cuprate superconductors. 
It is worth noting that in the literature the discussion concerning the occurrence of BKT physics in these layered 
anisotropic superconductors has been often associated to a somehow related issue, i.e. the nature of the pseudogap 
state above T c . Indeed, despite the intense experimental and theoretical research devoted to it, no consensus has been 
reached yet concerning its origin, with two main lines of interpretation based either on a preformed-pair scenario, or 
on the existence of a competing order, associated to fluctuations in the particle- hole channel 26 . The preformed-pair 
scenario has in turn triggered the attention on the role of SC phase fluctuations, which are expected to be very soft 
in these materials having a small superfluid density. Finally, the experimental observation of a large Nerst effect 33-35 
and diamagnetism 35 ~ 37 well above T c has been used to support the notion that phase fluctuations have a vortex-like 
character, as it is indeed the case within the BKT picture. However, the overall interpretation of the experimental 
data in terms of BKT physics is not so straightforward, despite several theoretical attempts based both on the 
mapping into the Coulomb-gas problem 38 ' 39 or on numerical simulations for the XY model. 40 On the other hand, 
a large Nerst effect arises also from the Fermi-surface reconstruction associated to stripe order 41 , or from ordinary 
GL fluctuations, 42 as observed for example in thin films of conventional superconductors 43 . Moreover, while the SC 
fluctuations contribution to the diamagnetism seems to fit the BKT behavior of the SC correlation length 36 , the 
paraconductivity shows usually more direct evidence of GL fluctuations 44 ' 45 . Also the issue of the superfluid-density 
jump is controversial: while it has been clearly identified in very thin films 13,16 ' 20 , it seems to be absent in bulk 
materials even for highly underdoped samples 46 , where the low superfluid density and the weak interplane coupling 
would make more plausible the presence of BKT physics. The aim of this paper is not to give a detailed overview of 
all the arguments in favor or against the occurrence of the BKT transition in cuprates, but to focus on the correct 
identification of the BKT signatures in a non-conventional quasi-2D superconductor. In general, a layered system 
with very low in-plane superfluid density and very weak interlayer coupling is one of the best candidates to observe 
those signatures of BKT physics that we mentioned above, i.e. a rapid downturn of the superfluid density coming 
from below T c and a regime of vortex-like excitations above T c . However, since the underlying superconductor is an 
unconventional one, the occurrence of the BKT physics could be masked by other effects, making its identification 
more subtle than in films of conventional superconductors. We notice that addressing this issue does not solve the 
more general problem of the nature of the pseudogap phase: for example, to estimate the extension in temperature 
of the regime of BKT fluctuations above T c one needs to consider physical ingredients that are beyond the BKT 
problem addressed here. Nonetheless, a deeper understanding of what could be the signatures of BKT physics can 
help discriminating its occurrence or not in unconventional superconductors. 

On the light of the above discussion, we will review in the present manuscript the work we have done in the last 
years to investigate the outcomes of the BKT transition in the presence of some additional ingredients that influence 
its occurrence in real systems without invalidating the basic physical picture behind it. The first ingredient that we 
will consider is the role played by the vortex-core energy fi, both in quasi-2D systems and in layered ones. As we shall 
see, while in the original formulation of BKT theory based on the XY model (1) /i is just a fixed constant times the 
coupling J, in real materials /i/ J can depend crucially on the microscopic nature of the underlying system. Since it 
represents the energy scale needed to create the core of the vortex, having "cheap" or "expensive" vortices can influence 
in a non-trivial way the tendency of vortex formation below and above the transition. Thus, while in 2D the critical 
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behavior will not change, in the layered 3D case the existence of expensive vortices can move the vortex-unbinding 
transition away from the temperature where it would occur in each (uncoupled) layer. The second aspect that we will 
discuss is the presence of inhomogeneity on a mesoscopic scale. This issue is somehow related to the effect of disorder 
on the BKT transition: however, instead of considering a model of microscopic disorder, we will implement a simpler 
approach where the spatial inhomogeneity of the superfluid density can be mapped in a probability distribution of 
the possible realizations of the superfluid-density values. This issue is in part motivated by several experimental 47-52 
and theoretical 53-55 suggestions that inhomogeneity on a mesoscopic scale can occur both in highly-disordered films 
of conventional superconductors 47-51 and in layered cuprate superconductors 52 , making then timely to investigate its 
effect on the BKT transition. Finally, we will discuss the role of a finite external magnetic field, an issue that is 
strictly related to the peculiarity of the BKT transition in a charged superfluid. In this case the motivation comes 
in part from recent experiments in cuprate superconductors, where anomalous non-linear magnetization effects have 
been reported 34-36 even for those samples 36 where apparently clear signatures of BKT physics appear. However, the 
main focus here is to establish a clear theoretical framework to deal with this complicated problem, based on the 
mapping into the sine-Gordon model. This last methodology aspect will serve as a general guideline for the present 
review paper. Indeed, we shall argue in the present manuscript that the mapping between the BKT transition in 
two dimensions and the quantum phase transition within the sine-Gordon model in one spatial dimension is the most 
powerful approach to explore the outcomes of the BKT physics beyond the standard results based on the original XY 
model. In particular, such a mapping provides us with a straightforward framework to explore the role of arbitrary 
values of fi/J, to include the effects of relevant perturbations (in the RG sense) as the inter layer coupling or the 
magnetic field, and to account at a basic level for the presence of inhomogeneities. Due to the several subtleties 
of such a mapping we shall first review the basic steps of its derivation in Sec. II, taking the point of view of the 
longitudinal vs transverse current decoupling in the XY model. Once established the general formalism and clarified 
the role of the vortex-core energy we shall address in Sec. Ill the consequences for the universal vs non-universal 
behavior of the superfluid density. In relation to the superfluid-density and paraconductivity behavior we shall give in 
Sec. IV a short account about the role of inhomogeneity and its observation in recent experiments in several systems. 
Finally, in Sec. V we give a detailed derivation of the sine-Gordon mapping in the presence of an external magnetic 
field, that completes our overview on the approach to the BKT physics in a real superconductor, and we discuss only 
briefly some physical outcomes related to the previous Sections. The concluding remarks are reported in Sec. VI. 



II. MAPPING ON THE SINE-GORDON MODEL AND THE VORTEX-CORE ENERGY 

As it is well known, the BKT transition occurs in three different physical phenomena, that belong to the same 
universality class: the vortex-unbinding transition within the XY model (1), the charge-unbinding transition in the 
2D Coulomb gas, and the quantum metal-insulator transition in the ID Luttinger liquid, as described by the sine- 
Gordon model. In the first two cases we are dealing with a classical model of point-like objects (vortices or charges) 
interacting via a logarithmic potential, that becomes short-ranged when the objects are free to move and to screen 
the interaction. In the latter case we deal with a quantum ID model, that becomes effectively a 2D one at T = 
where dynamic degrees of freedom provide the extra dimension. Even though all these analogies have been reviewed 
several times in the literature (see for example [8,9,56] just to mention a few references), nonetheless we will recall in 
this Section the main steps of the mapping between the classical XY model and the quantum ID sine-Gordon model. 
We shall take as a starting point of view the separation between longitudinal and transverse excitations of the phase, 
and we shall derive as an intermediate step the mapping on the Coulomb-gas problem, to make the physical aspects 
of the problem more evident. Instead in Sec. V, where we discuss the case of a finite magnetic field, we shall use a 
more formal approach, that has however the great advantage to provide us with an elegant and powerful formalism 
to discuss the case of a superconductor embedded in an external electromagnetic potential. 

As a starting point we shall consider a low-temperature limit for the XY model (1), where one could expect that 
the difference in angle between neighboring spins varies very slowly on the scale a of the lattice, so that one can 
approximate 9i — 6 i+ $ ss add(r)/d5 where 6(r) is a smooth function and 5 = x,y. Moreover, by retaining the leading 
powers in the phase differences from the cosine in Eq. (1) we find that in the low-temperature phase the model reduces 
to: 

Hxy = ^J dr(Ve) 2 = ±j rfrj 2 (r), (6) 

where we introduced, in analogy with the case of the superfluid, a current proportional to the phase gradient, j = V0. 
Because of the smoothness assumption the approximation (6) accounts only for the longitudinal component jy of the 
current, while vortices, i.e. singular configuration of the phase, can be associated to a transverse current component 
Indeed, a vortical configuration for the phase 9i of the XY model (1) corresponds to a non- vanishing circuitation 
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of the phase gradient along a closed line, that is nonzero only for a transverse current 

£j-d£ = y*(Vxj)-ds = J (Vxj 1 )-ds = 2^^ ft , (7) 

where = ±m is the vorticity of the i-th vortex, with m integer. We can then decompose in general the current of 
Eq. (6) as j — j|| + where V x j y =0 and V ■ jj_ = 0. One can easily see that the mixed terms J drj\\ ■ jj_ = in Eq. 

(6) vanish, so that longitudinal and transverse degrees of freedom decouple H = Hu + H±, and we can focus on the 
term H± = J drj^ to describe the interaction between vortices. By introducing a scalar function W the transverse 
current can be written as jj_ = V x (zW) = (d y W, —d x W, 0), so that V x jj_ = (0, 0, — V 2 W) and inserting it into Eq. 

(7) we conclude that W must satisfy the equation: 

\7 2 W(r) = -2np(r), p(r) = £ Qi 6(r - n). (8) 

i 

Eq. (8) is exactly the Poisson equation in 2D for the potential W generated by a distribution of point-like charges qt 
at the positions r^. Its solution is in general: 



W(r) = 2n J dr'V{r-r')p(r'), 



(9) 



where V(r) is the solution of the homogeneous equation for the Coulomb potential in 2D 



/r7k p ik r 
(2^1^' (10) 

so that V(r) ~ — lnr at large distances. Thanks to the results (8)-(9) H± can be written as: 

H± = J - J drj\ = ^ J dr{V x ZW) 2 = J - J dr(VW) 2 = ~J dvWV 2 W = 
= nJ J drW(r)p(r) = 2tt 2 J J drdr'p(r)V(r - r')p(r') = 

- 2n 2 jY / M 3 V{r l -r 3 ). (11) 

ij 

Eq. (11) expresses the electrostatic energy for a Coulomb gas with charge density p(r), completing thus the analogy 
between the system of vortices and the system of charges. The 2D Coulomb potential (10) shows the characteristic 
infrared divergence which reflects on the divergence of V(r = 0) in the thermodynamic limit, leading to the neutrality 
constraint for the gas. Indeed, by close inspection of Eq. (10) one sees that V(r = 0) = j 1 j L dk(l/k) <~ InL — > oo as 
L — >• oo. If we separate this divergent term by defining 

V(t) = V(0) + G{t), (12) 

where now G(r = 0) = 0, in Eq. (11) we obtain: 

2ir 2 J ^^QiQj[V{0) + G(vi — Vj)\ = 2tt 2 JV(0) fe*) + 2-7T 2 J ^ q l q J G{r l - Tj). (13) 
ij \ i / ij 

Since the Boltzmann weight of each configuration is e~P H± - , the divergence of V(0) in the thermodynamic limit leads 
to a vanishing contribution to the partition function, unless 

I>=0, (14) 

which means that only neutral configurations are allowed. A second consequence of the above discussion is that one 
should include a cut-off for the smallest possible distance between two vortices. Starting from the lattice XY model 
(1) a natural cut-off is provided by the lattice spacing a in the original model, which translates in the correlation length 
£ when applied to SC systems. The exact form of the function G(r) at short distances defines then the energetic 
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cost of a vortex on the smallest scale of the system, i.e. the so-called vortex-core energy. By computing G(r) on the 
lattice (see also Eq. (60) below), one can see that that at distances r > a it can be well approximated as 



G(r) ~ -\ -^-ln(-) . 
w 4 2f \aJ 



(15) 



Using the neutrality condition (14), the fact that G(0) = (so that in the last term of Eq. (13) one can use i ^ j) 
and the form (15), Eq. (11) can be written as: 



H ± = 27r 2 jJ2 q^G(r t -r J ) = -2n 2 jJ2 \ + i- In (^) 



where we used J2i^j QiQj = — Si 9? from Eq. (14) and we identified the vortex-core energy [i with 

7T 2 



- 2 J 



(16) 



(17) 



Finally, we can use the neutrality condition (14) by imposing that there are n pairs of vortices of opposite vorticity. 
Moreover, we shall consider in what follows only vortices of the lower vorticity Qi = ±1, so that reads: 



2n 

H± = 2nn — 7rJ^^ In ) £ i £ ji e» = ±1- 



(18) 



Eq. (18) describes the interaction between vortices in a given configuration with n vortex pairs. In the partition 
function of the system we must consider all the possible values of n, taking into account that interchanging the n 
vortices with same vorticity gives the same configuration (so one should divide by a factor l/(n!) 2 ). In conclusion Z 
reads: 



oo _^ 
n=l v ' 



dvi . . . dr 2n e 



oo „ 



where we introduced the vortex fugacity 



(19) 



(20) 



The explicit derivation of the partition function Z has the great advantage to allow us to recognize immediately 
the analogy with the quantum sine-Gordon model, defined by the Hamiltonian: 



Hs3 = ^l 



dx 



1 <~ .n2 9v 



K(d x ey + - (d x <t>y - ^cos{24>) 



(21) 



where 9 9 and 9 X represent two canonically conjugated variables for a ID chain of length L. with [0(x'), d x (p(x)] = 
m5{x' — x), K is the Luttinger-liquid (LL) parameter, v s the velocity of ID fermions, and g u is the strength of the 
sine-Gordon potential. In this formulation, the role of the SC phase is played by the field 9. Indeed, when the coupling 
g u = one can integrate out the dual field <f> to get the action 



K f 

S = — / dxdr [(d x 



9f + {d T ef 



(22) 



equivalent to the gradient expansion (6) of the model (1), once considered that the rescaled time r — > v s t plays the role 
of the second (classical) dimension. The dual field <f> describes instead the transverse vortex-like excitations. This can 
be easily understood by considering the quantum nature of the operators within the usual language of the sine-Gordon 
model. Indeed, a vortex configuration requires that § V9 = ±2ir over a closed loop, see Eq. (7) above. Since <fr is the 
dual field of the phase 9, a 2ir kink in the field 9 is generated by the operator e l2< ^, 9 i.e. by the sine-Gordon potential 
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in the Hamiltonian (21). More formally, one can show that the partition function of the <j> field in the sine-Gordon 
model corresponds to the (19) derived above. To see this, let us first of all integrate out the 9 field in Eq. (21), to 
obtain 

SsG= 2^kJ dr(V0)2 ~ f / drcos ( 2 ^' ( 23 ) 

The overall factor Zu = II g> o(l//3Jq 2 ) due to the integration of the 9 field (corresponding to the longitudinal 
excitations Zy = J T>6\\ e _/3i? H in Eq. (6) above) will be omitted in what follows. We can treat the first term of the 
above action as the free part So, and we can expand the exponential of the interacting part in series of powers, so 
that 

Z = [ V<\>eT Sa V -dv x ...dr p (-Y cos(20(ri)) • • ■ cos(20(r p )). (24) 
Here / T><j> is the functional integral over the 4> field. When we decompose each cosine term as 

i0( ri ) , p -i<t>(Ti) p i£0(ri) 

cos(20(r,)) = = J2 o ' (25) 

e=±l 

we recognize that in Eq. (24) we are left with the calculation of average value of exponential functions over the 
Gaussian weight So, i.e of factors 

( e 2i£ 4 e 4 *(r ( )) =^2^4^)^. (26) 

Here we used the well-known properties of Gaussian integrals 9 ' 56 that impose that the above expectation value is 
non zero only for neutral configurations Y7i=\ £ i = 0; m mu analogy with the result found above for the vortices. 
We then put again p — 2n. Taking for instance ei,...e n — +1 while e n +i, ■ ■ ■ £2n = —1 the combinatorial prefactor 
1 = p\ = l/(2n)! in Eq. (24) should be multiplied times the number ( 2 ^) = (2n)!/(n!) 2 of possibilities to choose the 
n positive £j values over the 2n ones. Thus, Eq. (24) reduces to: 



2 - Enp^r/*-*^^ < 27 » 



^ (n!)= 

1—1 K ' 

By comparing Eq. (19) and Eq. (27) we see that the vortex problem (as well as the Coulomb-gas problem) is fully 
mapped into the sine-Gordon model, provided that we identify: 

K=^ t (28) 
g = 27re-^. (29) 

As it is clear from the above derivation, within the XY model there exists a precise relation (17) between the value 
of the vortex-core energy /i and the value of the superfluid coupling J. This is somehow a natural consequence of the 
fact that the XY model (1) has only one coupling constant, J. Thus, when deriving the mapping on the continuum 
Coulomb-gas problem (18) p, is fixed by the short length-scale interaction, that fixes the behavior of G(r) in Eq. (15) 
and consequently the vortex-core energy value (17). In contrast, within the sine-Gordon language \i is determined 
by the value of the interaction g for the model (23), which can attain in principle arbitrary values. Thus, such a 
mapping is the more suitable one to investigate situations where fi actually deviates from the XF-model value, as it 
is suggested by the physics of various systems. 

A typical example is provided by the case of ordinary films of superconductors, where usually the BCS approx- 
imation -and its dirty-limit version- reproduces quite well the values of the SC quantities, like the gap and the 
transition temperature. 14,15 In this case, one would also expect that \x has a precise physical analogous with the loss 
in condensation energy within a vortex core of size of the order of the coherence length £ 0j 

H = n^£ cond (30) 

where e con d is the condensation-energy density. In the clean case Eq. (30) can be expressed in terms of J s by means 
of the BCS relations for e con d and £ - Indeed, since e con d — dN(0)A 2 /2, where N(0) is the density of states at the 
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Fermi level and A is the BCS gap, and £o = £,bcs — hvp/nA, where vp is the Fermi velocity, and assuming that 
n s = n at T = 0, where n — 2N(0)v F m/3, one has 

irh 2 n s d 3 3 . . 

^bcs = — : o = nJs—n - 0.95J S , (31) 

so that it is quite smaller than in the AY"-model case (17). This can have profound physical consequences on the 
manifestation of the BKT signatures in real materials, as we shall discuss in more details in the next Section. 



III. THE UNIVERSAL JUMP OF THE SUPERFLUID DENSITY 



The above derivation of the mapping between the XY model and the sine-Gordon model can be used directly 
to identify the two relevant running couplings K and g that must be considered under renormalization group. The 
coupling K (28) is connected to the superfiuid behavior of the system, thanks to the identification (5) of J with the 
superfluid stiffness of the system. Such an identification relies on the analogy between Eq. (6) and the kinetic energy 
of a superfluid. As a consequence, when also vortex excitations are present the physical J s must account also for 
vortex-antivortex pairs at short distances. 7 This physical picture has a precise correspondence on the values of the 
coupling constants under RG flow, whose well-known equations are 3,8,9 : 

% - < 2 - A ->»- < 33 » 

where I — lna/a is the rescaled length scale. The superfluid stiffness is then identified by the limiting value of K as 
one goes to large distances, i.e. 

TT 

Even though the behavior of the RG equations (32)- (33) has been described at length in several papers, we want to 
recall here the basic ingredients needed to describe the BKT transition. There are two main regimes: for K > 2 the 
r.h.s. of Eq. (33) is negative, so that g — > and K tends to a finite value — > K* that determines the physical stiffness 
J s , according to Eq. (34). Instead for K < 2 the vortex fugacity grows under RG flow, K in Eq. (32) scales to zero, 
and J s = 0. The BKT transition temperature is defined as the highest value of T such that K flows to a finite value. 
This occurs at the fixed point K = 2, g = 0, so that at the transition one always have 

K(l -+ 00, T BKT ) = 2, =► vJ fBKT) = 2 (35) 

J- BKT 

As soon as one goes to temperatures larger than Tbkt K — > 0, so also J s — > 0. As a result, one finds J(T^ KT ) = 0, 
i.e. the superfluid density jumps discontinuously to zero right above the transition. The equation (35) describes the 
so-called universal relation between the transition temperature Tbkt and the value of the superfluid stiffness J s at 
the transition, and represents a more refined version of the relation (3) based on the balance between the the energy 
and the entropy of a single- vortex configuration. 

It should be noticed that the BKT RG equations account only for the effect of vortex excitations, so that any 
other excitation that contributes to the depletion of the superfluid stiffness must be introduced by hand in the initial 
values of the running couplings. For example, in real superconductors there are also quasiparticle excitations, while 
in the XY model there are also longitudinal phase fluctuations, that give rise to a linear depletion to the superfluid 
stiffness Jo(T) = J(l — T/AJ). As a consequence the hand- waving argument usually adopted in the literature to 
estimate Tbkt in a system that is expected to have a BKT transition is to look for the intersection between the 
universal line 2T/7T and the J(T) expected from the remaining excitations except than vortices. However, such a 
procedure can only be approximate, since in the relation (35) the temperature dependence of J S (T) is determined also 
by the presence of bound vortex-antivortex pairs, which can renormalize J s already below Tbkt- This effect is usually 
negligible when /i is large, as it is the case for superfluid films 57 or within the standard XY model. However, as fi 
decreases the renormalization of J s due to bound vortex pairs increases, and consequently Tbkt is further reduced 
with respect to the mean-field critical temperature T c . As an example we show in Fig. 1 the behavior of J S (T) using 
a bare temperature dependence as in the XY model and switching the vortex-core energy from the value (17) to 
values smaller or larger. As one can see, for decreasing \i the effect of bound vortex-antivortex pairs below Tbkt is 
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FIG. 1: (a) Temperature dependence of the superfluid density for different values of the ratio fi/J, measured in units of the 
value (17) it has within the XY model. Notice that for small fi values the deviation from Jo(T) starts much before (blue arrow) 
than the temperature where the universal jump (red arrow) occurs. This behavior is indeed observed in thin NbN films, panel 
(b). Here we report data from Ref. [15] along with the theoretical BKT fit. Notice that the jump here is further smeared out 
by the inhomogeneity, see Sec. IV below. 



significantly larger, moving back the transition temperature to smaller values. The very same effect has been recently 
observed in thin films of NbN 14,15 , where it has been shown experimentally that the deviation of J S {T) from the 
BCS curve starts significantly below the transition temperature. Interestingly, the direct comparison between the 
experimental J s (T) and the results of the RG equations allowed the authors to show that the vortex-core energy in 
this system attains indeed a value of the order of the BCS estimate (31). 

It must be emphasized that the case of thin NbN films must be seen as a paradigmatic example of manifestation 
of the universal relation (35) , despite the fact that Tbkt is lower than expected for a standard view based on the 
AF-model results. Indeed, in this system once that the renormalized stiffness is of the order of 2T/ir the transition 
actually occurs. A different behavior is instead observed in the case we discussed recently 58 within the context of 
cuprate superconductors, that can be well modeled as weakly coupled 2D layers. In this case, an additional energy 
scale exists, i.e. the Josephson coupling J c between layers, which is also a relevant coupling under RG flow, so that a 
vortex-core energy different from the XY- model value can lead to a qualitatively different behavior. We notice that 
the case of weakly-coupled layered superconductors has been widely investigated in the past within the framework of 
a layered version of the XF-model (1). In this case the presence of a finite interlayer coupling J c « J cuts off the 
logarithmic divergence of the in-plane vortex potential at scales ~ a/^/yy, 59 where r] = J c /J, so that the superfluid 
phase persists above Tbkt, with T c at most few percent larger than Tbkt- 60 62 As far as the superfluid density is 
concerned, there is some theoretical 60 and numerical 63 evidence that even for moderate anisotropy the universal jump 
at Tbkt is replaced by a rapid downturn of p s {T) at a temperature T& larger than the Tbkt of each layer, but still 
very near to it T& ~ Tkt- Once more this result must be seen as an indication that the (only) scale J dominates the 
problem: the result can be different when fj, is allowed to vary, making the competition with the interlayer coupling 
more subtle. 

The analysis of the more general case has been done in a very convenient way in Ref. 58 within the framework of 
the sine-Gordon model (21), that must be suitable extended to include the interlayer coupling. As we said, in the 
sine-Gordon model (21) the variable 8 represents the SC phase. Since the phases in neighboring layers are coupled via 
a Josephson-coupling like interaction, the most natural assumption is an additional cosine term of strength J c for the 
interlayer phase difference, which translates in an interchain hopping term in the language of the ID quantum model. 
A similar model has been also derived recently in Ref. [64] by using as the starting point the Lawrence-Doniach model 
for the layered superconductor. The full Hamiltonian that we consider is: 58 

H = H S g[(f> m ,0 m ] — V / dscos[0 m -0 fn ,] > (36) 

m (m,m'> 

where m is the layer (chain) index and gj a = irJ c /T. In Ref. 58 we derived the perturbative RG equations for the 
couplings of the model (36) by means of the operator product expansion, in close analogy with the analysis of Ref. 
[65,66] for the multichain problem. Under RG flow an additional coupling g± between the phase in neighboring layers 
is generated: 





FIG. 2: (a) Temperature dependence of the superfluid stiffness J 3 (T) in the layered 3D case (symbols), taken from Ref. [58]. 
Here Jo(T) — J(l — T/4J) (solid line) mimics the effect of longitudinal excitations within the XY model. It is also shown for 
comparison the 2D case (lines), see Fig. 1. For /i < fixy Js(T) = at Td ~ Tbkt- As fi increases T<j increases as well, so that 
at Tbkt no signature is observed in J S (T) of the jump present in the 2D case, (b) Vortex-core energy in units of fixy as a 
function of the critical temperature of several YBa2Cu3C>6+:r samples, as derived in Ref. [67]. 



which contributes to the superfluid stiffness K s , defined as usual as the second-order derivative of the free energy with 
respect an infinitesimal twist S of the phase, d x 6 m — > d x 8 m — S. Thus, one immediately sees that Eq. (37) represents 
an interlayer current-current term, which contributes to the in-plane stiffness J s as: 

K s = K-nKg u J s = M±™E, (38) 

where n — 2 corresponds to the number of nearest-neighbors layers. The full set of RG equations for the couplings 
K,K s ,g u ,g Jc reads: 



dK 

~dl 

dg u 

de 

dK s 



= 2gj-K 2 gl (39) 

= (2-K)g u , (40) 

= -glKl (4i) 

9 J,- (42) 



dgj c _ ( 2 - — - Ks \ 
dl \ AK 4K 2 J 

Observe that for gj c — the first two equations reduce to the standard ones (28)- (29) of the BKT transition, and K s 
coincides with K. Instead, as an initial value gj c ^ is considered, the interlayer coupling increases under RG, 60-62 
contributing to stabilize the K parameter in Eq. (39), with a consequent slowing down of the increase of the g u 
coupling in Eq. (40). As in the pure 2D case, in the regime where K goes to zero g u increases, see Eq. (40), and 
vortices proliferate. However, in contrast to the single-layer case where g u becomes always relevant near K ~ 2, here 
thanks to the gj c term in Eq. (39), K can become smaller than 2 before than g u starts to increase. Since K s is 
controlled by the g u coupling alone, this means that the system remains superfluid in a range of temperature above 
Tbkt that depends on the competition between vortices (g u ) and interlayer coupling (g,j r ). 

The resulting temperature dependence of the superfluid stiffness J s for different values of the vortex-core energy 
is reported in Fig. 2a, where we added for the sake of completeness a temperature dependence of the bare coupling 
Jo(T) = J(l — T/4J), which mimics the effect of long-wavelength phase fluctuations in the XY model (1). Here 
the lines represent the pure 2D case (already shown in Fig. la above), while the symbols are the result of the RG 
equations (32)-(42) for a fixed value J c /J = 10~ 4 of the interlayer coupling. As one can see, as soon as a finite 
interlayer coupling is switched on, the jump of J S (T) at Tbkt disappears and it is replaced by a rapid bending of 
J S (T) at some temperature Td- However, while for \x < \xxy Td coincides essentially with Tkt, for a larger vortex-core 
energy Td rapidly increases and deviates significantly from the temperature scale where J a (T) intersects the universal 
line 2T/n. 

These results offer a possible interpretation of the measurements of superfluid density reported in strongly- 
underdoped samples of YBa 2 Cu306+2; 46 . These systems are exactly the ones where a BKT behavior could be expected: 
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they have a very low in-plane superfluid density and a large (77 ~ 10~ 4 ) anisotropy of the coupling constants. How- 
ever, the measured J S {T) goes smoothly across the Tbkt estimated from the 2D relation (35). In view of the above 
discussion, such a behavior does not automatically rules out the possibility that any BKT physics is at play: indeed, 
if the vortex-core energy is larger than in the XY model the transition can move away from the universal 2D case. 
Moreover, the effect of inhomogeneity can further round off the downturn induced by vortex proliferation (see next 
section), making it hardly visible in the experiments. In analogy with the case of conventional superconductors dis- 
cussed above, a hint on the realistic value of the vortex-core energy in cuprate systems can be inferred by the direct 
comparison between superfluid-density measurements and the theoretical prediction. We carried out such an analysis 
in the case of bilayer films of underdoped YBCO for several doping values, showing that \i attains always values 
larger than in the XY model, with 1.5 < fj,/ [ixy < 3, see Fig. 2b. Moreover, both [i and the inhomogeneity increase as 
the system gets underdoped, supporting further the interpretation that a similar effect can be at play also in the data 
of Ref. [46]. Interestingly, a very similar trend has been reported recently in NbN systems as disorder increases 15 , and 
it has been interpreted by the authors as an effect of the increasing separation between the energy scales associated to 
local pairing and superfluid stiffness as disorder increases. These results suggest that also disorder and inhomogeneity 
play a crucial role in the understanding of the BKT transition, as we shall discuss in the next section. 



IV. INHOMOGENEITY 



As we mentioned in the introduction, detailed measurements of superfluid density in thin films of superconductors 
became available only recently thanks to the efficient implementation of the two-coils mutual inductance technique, 
which gives access to the absolute value of the penetration depth of thin SC materials. Quite interestingly, in all 
the cases reported so far in the literature, concerning both conventional superconductors 14,15 and high-temperature 
superconductors 12 ' 13,16 ' 20,30 the BKT transition is never really sharp. At first sight one could wonder if any finite- 
size effect is at play, as due to several factors: (i) the existence of a finite screening length A ~ 2\ 2 /d due to the 
supercurrents in a charged superfluid 6 ; (ii) the finite dimension of the system or (iii) the finite length r u intrinsically 
associated to the probe, which uses an ac field at a typical frequency ui of the order of the KHz. In all the above 
cases one should cut-off the RG equations (32)-(33) at a finite scale £ max = \nL max /ao, leading to a rounding of the 
abrupt jump of the stiffness at Tbkt- However, in practice such rounding effects are hardly visible, since the decrease 
of K(£) at T > Tbkt is very fast, leading to visible rounding effect only for very short cut-off length scales of order of 
£ ~ 2 — 3. In real systems the cut-off length scales are usually much larger: for example, both in the case of the NbN 
films of Ref. [14] and in the case of thin YBCO films from Ref. [13] 1/A 2 is or order of about 1 /im -2 near the transition 
and d is of the order of 1 nm, so that the Pearl length A is of the order of 1mm, i.e. comparable to the system size, 
leading to £ max around 10 (a a is of the order of the coherence length £0 ~ 1 nm), which is practically an infinite 
cut-off for the RG. At the same time r u = y/lAD/ui is the maximum length probed by the oscillating field, where 
D ~ h/m = 10 16 A 2 /s is the diffusion constant of vortices and u) ~ 50 KHz is the frequency of the measurements 13,14 , 
giving again a large cut-off scale r w <~ 0.1 mm and no visible rounding effect, see Fig. 3a. It is worth noting that 
in the case of the experiments of Ref. 13 a second indication of the existence of pronounced rounding effects come 
from the experimental observation of a wide peak in the imaginary part of the conductivity around Tbkt- Following 
the dynamical BKT theory of Ambegaokar et al. 68 and Halperin and Nelson 5 such a peak is due to the bound- and 
free-vortex excitations contribution to the complex dielectric constant s(u>) = e 1 + is 2 which appears in the complex 
conductivity a = u\ + i(T 2 at a finite frequency cj: 

<r(u>) = . -. y (43) 

Here /io is the vacuum permettivity and we used MKS units as in Ref.s [13,67]. As in the case of the superfluid-density 
jump, we estimated 67 the width AT W in o\ (ui) due to the finite frequency uj of the measurements and we showed that 
it is expected to be much smaller than the experimental observations in Ref. 13, see Fig. 3a. 

On the light of the above discussion, we proposed in Ref. 67 that a more reasonable explanation for the rounding 
effects in the superfluid density and the large transient region in ai(iv) observed in the experiments is the sample 
inhomogeneity. In the case of cuprate systems such inhomogeneity is suggested by tunneling measurements in other 
families of cuprates, 52 where Gaussian-like fluctuations of the local gap value are observed. Recently similar local- 
gap distributions have been reported also in tunneling experiments in thin films of strongly-disordered conventional 
superconductors 47-51 , showing that an intrinsic tendency towards mesoscopic inhomogeneity appears even for systems 
with homogeneous intrinsic disorder. Even though the issue of the microscopic origin of this effect is beyond the scope 
of our work, nonetheless we found that a Gaussian-like distribution of the superfluid-stiffness Jo values around a given 
J can account very well for the data of Ref. [13] on two-unit-cell thick films of YBa2Cu306+ir- Thus, we compared 
with the experiments the quantity Ji n h(T) = J dJoP(Jo)J(T, Jo), where each J(T,Jq) curve is obtained from the 
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FIG. 3: Role of inhomogeneity, from Ref. [67]. (a) 1/A 2 and hqujoi evaluated at lo — 
parameter values appropriate for YBa2Cu306+a; films from Ref. [13]. Here p = 3(ixy- 



50 KHz for a single J(T) curve using 
The finite frequency leads to a sharp 



but continuous decrease of 1/A across Tbkt, along with a peak in eri. (b) 1/A and noujai evaluated at finite frequency taking 
into account the presence of inhomogeneities, as explained in the text. 



2D RG equations using a bare superfluid stiffness J = Jq — aT 2 . Each initial value Jo has a probability P(Jq) — 
exp[— (Jo — Jo) 2 /2cr 2 ]/(v / 27r(r) of being realized, where the bare average stiffness follows the typical temperature 
dependence due to quasiparticle excitations in a (disordered) d-wave superconductor, i.e. J(T) = Jo — aT 2 , with 
Jo and a fixed by the experimental data at low T, where J exp (T) is practically the same as J(T). Using a variance 
a = 0.05 Jo we obtained a very good agreement with the experiments near the transition, as far as both both the tail 
of A~ 2 and the position and width of <7i(w) are concerned, see Fig. 3b. 

Even though our approach to the issue of inhomogeneity is quite phenomenological, it has been shown more recently 
that it accounts very well also for the superfiuid-density behavior in films on NbN 14,15 . In all these cases indeed the 
rounding effect due to inhomogeneity is far more relevant than any finite-size effect due to screening or finite- frequency 
probes, explaining the lack of a very sharp jump even for relatively weakly-disordered films (see Fig. 1(b) above). 
Quite interestingly a similar physical picture seems to be relevant also for a completely different class of materials, 
i.e. the SC metal-oxides interfaces, where experimental data for the superfluid density are not yet available, but the 
nature of the SC transition can be investigated in an indirect way by means of the analysis of the paraconductivity 
above Tbkt 21 - For these systems as well a BKT transition is likely to be expected, since the SC interfaces are very 
thin 21 , of the order of 15 nm, which is much lower than the SC coherence length £o ~ 70 nm. Moreover, despite what 
occurs in thin films of conventional superconductors 15 , the drop of the resistivity above the transition is very smooth, 
suggesting that inhomogeneities on a mesoscopic scale broaden considerably the transition, as we discussed in details in 
Ref. [32]. On this respect, it is worth recalling that any analysis of the BKT transition based on the paraconductivity 
data alone can suffer of the unavoidable lack of knowledge about the exact extension of the BKT fluctuation regime. 
Indeed, as we mentioned in the introduction, the contribution of SC fluctuations to the paraconductivity is encoded in 
the temperature dependence of the correlation length £(T) both within the BKT and GL theory, see Eq. (4). Thus, if 
the transition has BKT character, one should expect a crossover from the BKT exponential temperature dependence 
of£(T), 

Zbkt = AexpQb/y/i), t= T ~ TbKT , (44) 

J- BKT 

to the power-law GL one, £gl ~ S,oT c /(T — T c ), where T c is the BCS mean-field critical temperature. In the case 
of thin films of superconductors the BKT regime is in practice restricted to a very small range of temperatures near 
Tbkt, due to the fact that T c at most twenty per cent larger than Tbkt- In the case of the SC interfaces it has 
been proposed instead in the recent literature 21,22,31 that T c is far larger than Tbkt, so that the whole fluctuation 
regime above Tbkt is dominated by BKT fluctuations alone, while the large tails observed experimentally near Tbkt 
should be ascribed to finite-size effects, see Fig. 4a. There is however a serious drawback in such interpretation, that 
originates from an incorrect application of the BKT relation (44) in the analysis of the experimental data. As we 
showed in details in Ref. [32] by means of a RG analysis of the correlation length, the parameter b which appears in 
the usual BKT formula (44) is directly connected to the distance t c — (T c — Tbkt) /Tbkt between the mean-field 
T c and the BKT transition temperature and to the value of the vortex-core energy with respect to the superfluid 
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stiffness: 



t 4 \l rr- . T C -T B KT 

Otheo ~ —o^j-Vtc, t c = . (45) 

7T Js 1 BKT 

As a consequence, b increases when the BKT fluctuation regime extends in temperature. In the typical fits of the 
resistivity proposed in Ref. 22,31 one obtains values of the b parameter that, according to Eq. (45) above, would imply 
a mean- field T c very near to Tbkt, see for instance Fig. 4a. In other words, the fit leads to values of b that contradict 
the a-priori assumption that the whole fluctuations regime is dominated by BKT fluctuation, as described by Eq. 
(44). Moreover, also the interpretation of the tails extending below Tbkt as due to finite-size effects 31 leads to 
unphysical low values for the size of the homogeneous domains. 32 One would have thus to invoque a BKT transition 
of a completely differente nature (dislocations of a vortex crystal 69 ), as it has been suggested in Ref. [21]. However, 
even in this case there is no theoretical understanding how to reconcile such contradictory numbers obtained in the 
analysis of the resi 




FIG. 4: Analysis of the paraconductivity effect on the resistivity of SC interfaces at the Tbkt transition. The experimental 
data for the resistivity R normalized to the normal-state value Rn are taken from Ref. [22]. Panel (a) shows an example of the 
approach proposed in [21,22,31], where the whole range of temperatures above Tbkt is dominated by SC fluctuations having 
BKT character, and finite-size effects are responsible for the resistive tail below Tbkt- However, the b value obtained by the 
BKT fit based on Eq. (44) for the correlation length implies, according to Eq. (45), the T c value marked by the blue arrow, 
invalidating thus the assumption itself that the whole fluctuation regime has BKT character. Panel (b) is taken from Ref. 70 
(Fig. 5.1 in that paper) and elucidates the approach developed in Refs. 32,70, based on the interpolation scheme (46) between 
BKT and GL fluctuation. The curve labeled "GL+BKT+inhom" has been computed by solving a random-resistor network 
problem with a Gaussian distribution (a = 0.035 K) of the critical temperatures centered around Tbkt = 0.2 K, marked by an 
arrow along with the corresponding T c — 0.22 K. The paraconductivity of each resistor follows Eq. (46) with b = 0.25, which 
according to Eq. (45) corresponds to t c = 0.1 when fi/ J s ~ 2, an intermediate value between the XY model (17) and BCS (31) 
estimates. The curve labeled "GL+BKT" corresponds to the homogeneous case with a transition temperature Tbkt- 



An alternative, and less far fetched, interpretation of the resistivity data can be based instead on the use of the 
interpolation formula for £(T) proposed long ago by Halperin and Nelson 5 

w £(T)=£oAsinhA. (46) 

Rn i + (£/£o) vt 

One can easily see that Eq. (46) reduces to £bkt for small t and to £cl for large t. Using the estimate (45) for the b 
parameter one realizes that the crossover occurs around t ~ b 2 oc t c . so that for (realistic) small values of t c most of the 
fluctuation regime is dominated by GL-type fluctuations. Moreover, we attributed the broadening of the transition to 
the effect of inhomogeneity, that induces a distribution of possible realizations of local R/ Rn values corresponding to 
the local J values discussed above. An example is shown in Fig. 4b, where the average resistivity has been computed 
by solving a random-resistor-network model for a set of resistors undergoing a metal-superconductor transition, as 
discussed recently in Ref. [70]. It must be emphasized that the issue of the physical value of the b parameter in the 
BKT expression (44) of the correlation length has been often overlooked in the literature also in contexts of other 
materials 24,28 , leading to questionable conclusions concerning the existence of a BKT physics based only on erroneous 
fits of the resistivity. In contrast, a recent analysis on NbN films 15 , based on the comparison between the estimate of 
b deduced from the analysis of the superfluid density below Tbkt and the resistivity above Tbkt has confirmed the 
theoretical estimate (45) of the b parameter, that must be used as an unavoidable constraint in any BKT fit of the 
paraconductivity. 
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V. THE CASE OF A FINITE MAGNETIC FIELD 



As we discussed in the previous Sections, the sine-Gordon model provided us with the most appropriate formalism 
to investigate not only the role of a vortex-core energy value different from the XY-mode\ one, but also the occurrence 
of BKT physics in the presence of a relevant perturbation (in the RG sense), as the interlayer coupling. It is then 
worth asking the question of the possible relevance of such a mapping for an other typical relevant perturbation, i.e. 
a finite external magnetic field. The nature of the BKT physics in the presence of a magnetic field has been of course 
already investigated in the past 4 ' 8 ' 71 , with a renewed interest in the more recent literature 38-40 triggered once more 
by experiments in cuprate superconductors carried out at finite field, as the measurement of the Nerst effect or the 
magnetization. Unfortunately, contrary to the case of the B = transition, the efforts have been partly unsatisfactory. 
In particular most of the literature on the subject rested on the use of the mapping into the Coulomb-gas problem, 
where the effects of the magnetic field can be incorporated as an excess of positive charges, in analogy with the finite 
population of vortices with a given vorticity due to the presence of the external field. There is however a fundamental 
drawback of this approach, since it gives the physical observables as a function of the magnetic induction B instead of 
the magnetic field H. This is of course not convenient at low applied field, since inside the superconductor B vanishes 
even for a finite external field H. Motivated by this observation and by the occurrence of an anomalous non-linear 
regime for the magnetization measured in cuprate systems, we showed recently 72 that a suitable mapping into the 
sine-Gordon model provides a very simple and physically transparent way to deal with the finite magnetic-field case, 
leading to a straightforward definition of the physical observables and clarifying the role of both B and H. Since the 
construction of the mapping has not been shown in Ref. [72] due to space limitation, we shall discuss here both the 
model derivation and its basic physical consequences. This will complete our overview of the sine-Gordon approach 
to the BKT physics, and it will allow us also to clarify some aspects of the charged superfluid that have not been 
underlined yet in the previous Sections. 

As a starting point we use again the XY model (1), where the coupling to B is introduced via the minimal-coupling 
prescription for the vector potential A, 

H = J V [1 - cos(0< - 9 3 - F i3 )} (47) 



<i,j> 



where 



2tt f r+IJ - 2ira 
Fij = F,(r) = — A ■ dl « -3-A„(r). (48) 

^0 Jr 

Here ji = x, y is the vector from one site to the neighboring one and $o = hc/2e is the flux quantum. As a first step 
we make use of the Villain approximation, that amounts to replace the cosine in Eq. (47) with a function having the 
same minima for each multiple 2inn of the gauge-invariant phase difference, i.e. 



m— oo 



exp[— /?J(1 — cos(6>i - 6j - Fij)} => ^ exp 

m— — oo 



2 



- 2-Kmf 



(49) 



By following the analogous derivation given in Ref. [56] for the case A = 0, we make use of the Poisson summation 
formula 



2-Kilz 



h(m) = / dzh(z)e 

— r>c> 1 = — r*~i 



(50) 



m— — oo 



and by performing explicitly the integration over z for each link ij we are left with the following structure of the 
partition function: 

/■27T ,2. 

Z = ZjJ2 d6 1 ...de N e~ T '<^>^ +iT '^i) lis{ei ~ ei ~ Fii \ (51) 
ih 3 } Jo 

where N is the number of lattice sites and Uj are 2N integer variables defined for each link ij. The prefactor 
Zj = [(1//?J) 1 / 2 ] 2A ' = (1//3J) N accounts for the z integration above for each link. Since each 9 variable is defined 
on a period, we have that d6e l0a = 5(a). As a consequence, the integration over the 6i variables in the above 
equation leads to N constraint equations: 

^(r)-^(r- M ) = 0, (52) 
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which are the discrete equivalent of V • 1 = 0. These equations can be satisfied by defining for each site of the lattice 
a field n such that: 

l x (r)=n(r)-n(r-y), (53) 
l y (r) = n(r — x) — n(r), (54) 

i.e. the discrete equivalent of the relation 1 = V x (nz). By means of Eqs. (53)-(54) Eq. (51) can be then rewritten 
as: 

Z = Zj e^,w^»M) 2 e - i ^S r -(A X An) ! (55) 

Mr)} 

where A M = n(r + /i) — n(r) is the discrete derivative. Finally we can use again the Poisson summation formula to 
get: 

Z = Zj[V(j){v) ^ 2h( A ^)) 2 +^ i T, r 4'(r)rn(r)-t^-E r z-(AxA4>)^ ^ 

{m(r)} 

where the m(r) variables assume arbitrary positive and negative integer values. Before completing the mapping into 
the sine-Gordon model we notice that in the above equation the (j> variables can be integrated out exactly. By rewriting 
the last term in the action as ^ r z ■ (A x A(j>) = J2 r <f>z • (A x A) = Ba J2 r § one can easily see that the final result 
is given by 

z = Z\\ er 2n20J ^y [m{r) - f]u< - r - r ' )[m( - r ' ) - f \ / = — (57) 

Mr)} $ ° 

Here Z^ = ZjH q>0 (f3J/c{ 2 ) = U q>a (l/ (3Jq 2 ) is the overall contribution due to the longitudinal excitations. In full 
analogy with the result of Sec. II, the longitudinal modes decouple from the transverse ones, so their contribution Z\\ 
to the partition function can be discarded in what follows. The function V(r) = ^2 k e t]i ' r U(k) is defined trough the 
Fourier transform J7" 1 (k) = (4 — 2cosfc. E — 2cosfcj,) of the A M operator on the square lattice, i.e. 

/d 2 k e ikr 
(2tt) 2 [4- 2 cos k x - 2 cosfc,,]' ^ 

Eq. (57) generalizes the Coulomb-gas formula (11) above to the case of a finite magnetic field. As we did in Sec. (II), 
we can separate in V(r) the singular part in r = by defining the regular function G(r), i.e. V(r) = V(0) + G(r). 
Then one sees that 

J2[m(r)-f]V(r-r')[m(r')-f} = 



V(0) 



-1 2 



+ Y^[m(r)-f}G(r-r')[m(r')-f] 



It then follows that also in this case only neutral configurations have a statistical weight different from zero. However, 
in the presence of a magnetic field the neutrality condition reads: 

J2(m(r)-f) = N v -—— = (59) 

which means that the total flux N v &o carried out by the (unbalanced) vortices equals the total flux Ba 2 N of the 
magnetic field across the sample. The definition (58) allows us also to determine the value (17) of the chemical 
potential \i in the lattice XY model. Indeed, one can see that at the scale of the lattice spacing V(r) — V(0) gives: 



d 2 k cosfc T 



) 2 [4 — 2 cos k x — 2 cos k y ] 

cos k x + cos k y — 2 1 
2 J (2tt) 2 [4 -2 cos k x - 2 cos k y ] ~~ ~4' ' ' 
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so that from Eq. (57) it follows that the cost to put two vortices at distance a apart is (3[i = (3ir 2 J/2, consistent with 
the value (17) that we quoted above. 

Let us now go back to the Eq. (56) and let us complete the mapping into the sine-Gordon model. We notice that 
in Eq. (56) the variable (f> is still defined on the square lattice. We can however resort to a continuum approximation 
by taking into account the energetic cost fi of the vortex creation on the shortest length scale of the problem via a 
chemical-potential like term e lnv ^ m ( r ) in Eq. (56): 

Z = [ V4>{r) ^ e~ Er ^ w( A ^( r » 2 + 2 «E r E r i-{Ax A« e l n!) m J (r) ^ 

{m(r)} 

The vortex-core energy term favors the formation of vortices of smallest vorticity, i.e. m(r) = 0, ±1. If we limit 
ourselves to this case the sum over the integer variables m(r) can be performed explicitly as: 

e in Vm 2 +2^ mi p = 1 + 2ycos(2ncf)) w e 2vcos{2lut,) (62) 

m(r)=0,±l 

Inserting this into Eq. (56), taking the limit of the continuum (^ r — > (1/a 2 ) J d 2 r), and rescaling <fi —> 4>/ 7T we finally 
obtain a partition function expressed in terms of the <j> field only, that generalizes Eq. (23) above: 



drdz 



(V</f>) 2 g nA 2% 



2ttK na 2 $ 



cos 20+ — A-(V x z<j>) 



o 



S(z), (63) 



where we used again the definitions (28)-(29) for K, g. In Eq. (63) we added also the explicit z dependence of the 
action, which is needed since the A field depends in general also on the z out-of-plane coordinate. The 8{z) function 
gives the proper boundary conditions for a truly 2D case (where there is no SC current outside the plane), while in the 
physical case of a SC film of finite thickness d we shall assume that the sample quantities are averaged over \z\ < d/2. 

The action (63) and its corresponding partition function Z = J V<f)e~ SB allow for a straightforward definition of the 
physical observables at finite field. For example, the electrical current follows as usual from the functional derivative 
of F = — (1//?) InZ with respect to the gauge field A, i.e 

T / \ OF 2ickBT ,_ .. . . A . 

Js(r ' z) = ~ c dM^) = — $^ <v x WWW- ^ 

The Eq. (64) makes it evident that the current is purely transverse, as expected for vortex excitations, according to the 
discussion given in Sec. II above. A second quantity that can be easily obtained is the magnetization M = (B — H)/4ir, 
defined as the functional derivative of F with respect to B(r, z) = V x A. By integrating by part and using the identity 
J (Ax V(f>) ■ z = J V</> ■ (z x A) = — J (fN ■ (z x A) = J 4>z ■ (V x A), we can rewrite the last term of Eq. (63) as 



/ 



drdz 



W) 2 9 _ OJ , , 2* 
'o 



2*K ^-^ + -B(r,,).^(r) 



6(z), (65) 



As a consequence the functional derivative with respect to B gives immediately 

and analogously the uniform magnetic susceptibility \ — 9M z /dB z is: 

X = -{Ak B T)/{d<g) ( dr [<0(r)0(O)) - (<^(r)) (0(0))] . (67) 



Notice that the total magnetic moment M. associated to the current J s is A4 = (l/2c) J drdz (r x J s ) 73 . Using the 
definition (64) of J s , one can easily verify that M = J d 2 rdzM., i.e. the magnetization is the density of magnetic 
moment of the sample, as expected 73 . Finally, by exploiting the fact that e _/3/i e ±l2< ^ is the operator which creates up 
and down vortices with density n± respectively, we have a straightforward definition of the average vortex number 
rip = a 2 ((n + ) + (n_)) and of the excess vortex number n = a 2 ((n + ) — (n_)) per unit cell as a function of (f> as: 

n F = 2e-^(cos(20)), n = 2e" /3 ' t (sin(20)). (68) 

In Eq. (66) above the average value of <j) is computed with the action (63), so that the magnetization is given 
as a function of the magnetic induction B. However, as we mentioned above, at low external field it would be 
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more convenient to compute M as a function of the applied field H. This can be achieved by using the Gibbs free 
energy Q = — k B T\nZ G , where the partition function Z G includes also the contribution of the electromagnetic field 
Zq = J V<j)VAe- s and: 



f f 
/ dvdz < - 



SB + ldr dz{ ^l A S- iV A X ^ U \- (69) 
1 87rfc B T 47rfc B T ' 



Before discussing explicitly the case of a finite external field H we would like to stress how Eq. (69) gives also a 
very convenient description of the role of charged supercurrents in a 2D superconductor. Indeed, even when H = 
the electromagnetic field A in Eq. (69) above describes the magnetic field created by the current themselves in a 
charged superfluid. In this case, since the SC currents live in the plane, also A is a two-dimensional vector, so that 
V x A = (—d z Ay,d z A x ,d x A y — d y A x ). Moreover, is we choose the Coulomb (or radial) gauge V • A = 0, we have 
that in Fourier space (k|| x A) 2 = fc 2 A 2 . We can then rewrite in Fourier space the terms in A of Eq. (69) at H = 
as: 



/ 



d 3 k { 2 (^+k 2 ) 1 



(2tt) 3 I $ " 11 "' 8nT 
By integrating out A at Gaussian level we obtain a (f> 2 contribution to the action of the form: 

d 3 k 8irT k 2 



w * §( « + k f)'*< k ''>i- ™ 

Since <\> depends on ky only, we can integrate out k z and obtain that the overall Gaussian action for the <j> field reads: 

rf 2 k„ 1 

where we defined: 



* G= / (2#2^ [fc " +fc " A_1]|(/)(k '' )|2 ' (?2) 



1 _ 8tt 2 KT _ d 
A ~ $ 2 = 2A 2 ' 



(73) 



The last equality follows from the definitions (28) and (5) of K = ttJ/T and J = $Q(i/167r 3 A 2 , respectively, and it 
allows us to identify A with the so-called Pearl screening length 6 . Indeed, due to the fc|| A -1 terms in Eq. (72), one sees 
that the potential between vortices, which according to the above discussion is the Fourier transform of the Gaussian 
(j) propagator, decays as e~ r ' A at scales r » A, instead of the usual log r dependence observed at all length scales in 
neutral superfluids. 

Let us go back now to the case of a finite external field in Eq. (69) and let us integrate again the gauge field A. 
Since now it is present an additional term — («/47rT)A(k) • (k x H(— k)) in S, one obtains the following contribution 
to the action: 



8ttT 



^(4xk||) + ^(kxH(-k)) 



(74) 



The quadratic term in <j> corresponds to Eq. (71) above, leading to the screening of the vortex potential. The remaining 
terms can be written as: 

d 3 k 2i 



and 



Using the identity 



/d^k 2? 
(27r)3 g ofe2 (* >< k ll<K k ll)) • (k x H(-k)), (75) 

d 3 k 1 (k x H) 2 



~>H-H 



(2tt) 3 8ttT k 2 



(76) 



d 3 k 1 r 1 

— Fl (k)F 2 (-k)- = J rf 3 rrf 3 r , Fi(r)jF2(r0 ___ 5 (77) 



(2 
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and the Maxwell equation relating the magnetic field H to the distribution of the external current 3 ext producing the 
field itself 

47T 

V x H = — 3 ext , (78) 

c 

one can easily see that Eq.s (75) and Eq. (76) can be written in real space as: 

J 47r|r — r'| 

and 



1 A rf 3 rd 3 r , J e ^ W'Jext(r') 

SttT 7 4vr|r-r'| ' v ' 

respectively. By integration by part Eq. (79) can be rewritten as: 

= A f rf 3 rrf 3 r ,^V^)].J rat (r') = 
c J |r — r 'I 

c J r - r'| 



d»rrfV*(r)* • [J - t( ° X( ^ r0 

r — r' P 



= |^|rf 3 r^(r)z.H (r). (81) 

In the last equality of Eq. (81) we introduced the reference field H°, which corresponds to the magnetic field generated 
by the same distribution of currents 3 ext in the vacuum. According to the Laplace formula, H° is given exactly by 

■w = \ i ^ j -y. x j|,-' ) . <m 



H' 

leading to Eq. (81) above. One can also recognize in the term (80) the magnetic energy density associated to the 
reference field H°, 



l - J d 3 r(H ) 2 . (83) 



S = 



DH ~ H= 8nT 

Indeed, since H° is the field created by the currents 3 ext in the vacuum it satisfies V • H° = 0, so that (k x H ) 2 = 
(H°) 2 fc 2 . Thus H°(k) 2 = (V x H°)/fc 2 = 3 ext (k) 2 /k 2 that is the Fourier transform of Eq. (80). In summary, the full 
sine-Gordon action after integration of the gauge field can be rewritten as: 

Eq. (84) is the desired result to be used to evaluate the physical observable as a function of the reference field H°. Once 
more, the sine-Gordon mapping turns out to provide a quite powerful framework for the investigation of the BKT 
physics of a superconductor embedded in an external field. Indeed, apart from the fact that it includes automatically 
the screening effect of the supercurrents discussed above, the action (84) expressed in terms of H° has two main 
advantages. First of all, H° is the field quoted in the experimental measurements, since what is known a priori 
are only the generating currents 3 ex t- Indeed, H° does not coincide in general with the real field H even outside 
the sample, since the H configuration takes into account also the field exclusion from the SC sample, the so-called 
demagnetization effects. For simple sample geometries one can include these effects in a demagnetization coefficient 
r], and write in general the following relation between B = H + 47rM and H 073 : 

(l-r ? )H + 77B = H°^B = H° + 47r(l-r7)M. (85) 
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In the complete Meissner phase one has B = 0, which implies — 47rM = H. However, from Eq. (85) it follows that 
H = H°/(l - n) so that: 

1 H° 

M = --- . (86) 

47T 1 — 7] 

While for a cylinder 77 = and H = H°, for a film of thickness d and transversal dimension R one has that 77 ~ 1 — d/R, 
so one expects to find M <~ (R/d)H.° below H c \, i.e. a much smaller critical field with respect to the same system in 
the 3D geometry 73 ' 74 . Since the magnetization M calculated from Eq. (84) is already a function of H°, it will include 
automatically all the demagnetization effects and the complications of the thin-film geometry. 

These properties have been derived in Ref. [72], where the magnetization has been computed by means of a 
variational approximation for the cosine term in the model (84). While we refer the reader to Ref. [72] for more 
details concerning these calculations, we would like to mention here one particular result, that is related to the 
discussion of the previous Sections. It concerns the behavior of the field-induced magnetization above Tbkt, that is 
expected 5 to be proportional to H with a coefficient depending on the SC correlation length: 

M = -*-§/"■ (87) 



In full analogy with the paraconductivity discussed in the previous Sections, the functional dependence of the low-field 
magnetization M on the BKT correlation length £ in Eq. (87) is the same as in the GL theory. While this result was 
already known in the literature 5 ' 38 , our calculations based on the model (84) allowed to establish an upper limit Hi 
for the validity of the linear regime (87) 



H<Hi = 0.1^^pSZ. (88) 

Notice that the above relation can be approximately expressed as the condition £ » £b for the low- field limit to 
be applied, where £ 2 B = <&q/H is the magnetic length scale. 43 As T approaches Tbkt £ increases rapidly and the 
field Hi becomes rapidly smaller than the lowest field accessible in the standard experimental set-up. This effect 
can explain for example the non-linear magnetization effects reported recently in several measurements in cupratc 
superconductors 34-36 . Indeed, the persistence of a non-linear magnetization up to H <~ 0.01 T in a wide range of 
temperatures above Tbkt can be a signature of the rapid decrease of Hi as T — > Tbkt, which does not contradict 
but eventually support the BKT nature of the SC fluctuations in these systems. Moreover, since £ increases as fi 
increases, the extremely low values of Hi measured in Ref. 36 suggest a value of fi larger than (j,xy, in agreement with 
the result discussed in Sec. Ill based on the analysis of the superfluid density. On the other hand also the existence of 
inhomogeneities can alter the straightforward manifestation of a linear magnetization above Tbkt, an issue that has 
not been explored yet neither in the context of cuprates nor in the case of conventional superconductors. Finally, we 
would like to mention that even though some theoretical work exists 38 on the RG approach to the BKT transition at 
finite magnetic field based on the Coulomb-Gas analogy, a full analysis of the more general model (84) is still lacking. 
Such an approach could eventually improve the estimate (88) of the linear regime, based on a variational calculation 
that is not expected to capture the correct critical behavior as the transition is approached. 



VI. CONCLUSIONS 



It is clear that BKT theory has profoundly changed our understanding of quasi-2D superconductors and given us a 
tool to tackle such challenging and interesting problems. However, more than 40 years after the original discovery, the 
occurrence of the BKT transition in several quasi-2D superconducting materials remains partly controversial. One can 
in general identify two possible sources of discrepancies between theoretical predictions and the current experimental 
scenario. From one side, the original formulation was based on the paradigmatic case of the XY model, that is 
only one possible model where the BKT transition occurs. Even though it correctly reproduces the critical behavior 
of all the systems belonging to the same universality class, quantitative discrepancies away from criticality can be 
observed in different models. This is the case of the strong superfluid-stiffness renormalization below the transition 
temperature Tbkt in the case of superconducting films of conventional superconductors, where the vortex-core energy 
attains values significantly different from the XY model prediction. From the other side, emerging new materials and 
improved experimental techniques offer new scenarios for the occurrence of the BKT transition, which coexists with 
several other phenomena. An example is provided by the case of cuprate superconductors, that are layered systems 
formed by strongly-correlated 2D SC layers. In this case, the deviations of the vortex-core energy from the XF-model 
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value can eventually lead to a qualitative different behavior of the superfluid-density jump at the transition or to strong 
non-linear field- induced magnetization effects above Tbkt- I n the present article we reviewed a possible approach to 
all these issues based on the sine-Gordon model. Even though this is certainly not a new approach for the pure 2D 
case, in the presence of additional relevant perturbations it provides a very convenient framework to investigate the 
BKT physics. Indeed, it allows not only to incorporate easily the effects of a vortex-core energy value different from the 
XY model, but also to describe the coupling to the electromagnetic field in a clear way, giving a straightforward and 
elegant description of the charged superfluid. Finally, we would like to emphasize once more that a quite interesting 
issue, that applies equally well to conventional and unconventional superconductors, is posed by the role of the intrinsic 
sample inhomogeneity. Even though we outlined here a kind of mesoscopic approach to the emergence of spatially 
inhomogenous SC properties, a more microscopic approach to the effect of disorder on the BKT transition would be 
required, as suggested by some recent numerical works 75,76 . The theoretical and experimental investigation of this 
issue will certainly offer an other perspective on the BKT transition in low-dimensional superconductors. 
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